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We discuss a possible origin of Tsallis' statistics from the correlation among constituents which 
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reduces the phase space of the system. We also show that in the system of coupled linear harmonic 
O ■ oscillators can exhibit a Tsallis type behavior. 
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I. INTRODUCTION 



If we want to determine precisely the dynamical behavior of a macroscopic system out 
of equilibrium, we need in general to deal with highly coupled equations of motion related 
to the full microscopic degrees of freedom of the system. However, when we are interested 
in knowing only the dynamics of macroscopic observables, we can sometimes expect a large 
amount of reduction of degrees of freedom relevant to the dynamics of these macroscopic 
observables. They should be described by a closed set of simpler coarse-grained equations. 
Such a situation happens when the space-time scales for our macroscopic observables are 
clearly separated from those of the microscopic ones. The hydrodynamics is such an example. 

Hydrodynamics is derived by assuming the local thermal equilibrium: at any point of a 
fluid, there exists a finite domain (fluid cell) in which the fluid can be regarded as homoge- 
neous and in thermal equilibrium, so that the usual thermodynamic relations hold for this 
piece of fluid. The physical properties of the fluid are determined by these thermodynamic 
relations. Furthermore, although finite, the size of the domain is considered to be negligible 
compared to the typical scale of the system, so that the thermodynamic variables can be 
considered as continuously varying functions of position and time. 

On the other hand, strictly speaking, the theromodynamic description of a matter is 
only possible when the time scale in question is much larger than the relaxation time of 
the matter needed to attain thermodynamic equilibrium. In other words, for hydrodynamic 
description, the fluid cell should have the minimum scale to realize local thermal equilibrium. 
For example, for a gas, such a domain should at least be several times larger than the mean- 
free path of the gas particles, otherwise fluctuations are too large and the thermodynamic 
description of the gas property will fail. It is thus reasonable to consider that hydrodynamic 
approach is not applicable when the size of the system is relatively small so that the typical 
space-time inhomogeneity scales are comparable or smaller than those required for the local 
thermalization of the system. In fact, it is believed that the deviations from the usual 
thermodynamics is an important factor in physics associated with molecular scales such as 
biological systems [lj]. 

A parameter commonly used to evaluate quantitatively the applicability of hydrodynamic 
approach is the Knudsen number K, defined as the ratio, K = X meso / Xhydro where A meso is 
the mesoscopic scale of the microscopic structure (in the case of a gas, it is the mean-free 
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path of the gas particles ) and Xhydro is the typical hydrodynamical inhomogeneity scale. If 
K <C 1 everywhere, a hydrodynamic description will be applicable as the dynamics of the 
system. This also means that the time scale of the macroscopic (hydrodynamic) scales are 
much larger than the thermodynamic relaxation time. In general, the hydrodynamic time 
scale is of the order of \h y dro/v soun d , where v soun d is the velocity of the sound of the fluid. 

On the other hand, there are several counter examples where hydrodynamic models seem 
to work for K values apparently not so small. One of these is the physics of relativistic 
heavy- ion collisions j^]. To analyse the behaviors of huge amount of produced particles in 
the relativistic heavy-ion collisions, we simply assume that the many-body dynamics can be 
approximately given by relativistic hydrodynamics (in particular, perfect fluid). Since the 
typical dimension of the whole system is initially from few fm (10 _15 m) (or even less if we 
consider the Lorentz contraction due to the relativistic effects) and expands violently. For a 
hydrodynamic description to be valid from the beginning of this expansion, we expect that 
the fluid cell size should at least be less than one order of magnitude of lfm. That is, we 
should use the cell size of the order of 0.1 fm, which is already comparable, if not smaller, to 
a typical value of mean-free path of the constituent particles, partons. The value, 0.1 fm is 
even closer to the scale of quantum correlation length. In spite of these facts, the collective 
motion of the matter produced in the relativistic heavy ion collision is well described by the 
hydrodynamic model. The hydrodynamic model is a very important tool to understand the 
experimental data. 

Then how should we understand such a situation? Here, we have to recall that the local 
thermodynamic equilibrium is a sufficient condition for the thermodynamic relations to be 
valid, but not a necessary condition. For example, let us consider an ideal Boltzmann gas. If 
the system is isotropic in momentum space at any position, the energy density and pressure 
are given as 

e(r) = j d 3 P E(p)f(p,r), (1) 

and 

P(r) = \ [ d 3 p (?■ if) f(p,r), (2) 
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where / (p,r) is the one-particle phase space distribution function and v = p/E. If these 



particles are massless (or ultra- relativistic), then \p\ = E (p) , so that we always have [11] 



e = 3P, (3) 
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independent of the form of the distribution function / 12]. 

As we saw above, what we need to arrive Eq.Q is only the condition of local isotropization 
of the particles in the momentum space, and not the thermal equilibrium. It may well be 
possible that the process of local isotropization in momentum space can be attained much 
faster than the complete thermalization. In fact, for a quark-gluon plasma which is now 
believed to be produced by relativistic heavy-ion collisions, such a dynamical mechanism 
has been studied in analogy with the plasma instabilities {3)]. This is an example where 
the thermodynamic relations become valid very quickly even the system is not in a real 
thermal equilibrium. It is also shown that the nonequilibrium dynamics of quantum scalar 
field exhibits such a prethermalization behavior J4]. These considerations indicate that there 
are situations where hydrodynamical description can still be valid even for a system where 
the usual criteria in terms of Knudsen number is apparently violated. From now on, we 
refer such a state of matter to "prethermalized" state, for which functional relations among 
macroscopic observables (such as equation of state) emerge, although the true thermalization 
is not yet attained. 

Now, two basic questions arise here. One is in which situations such a prethermalized 
state occur, and the second, whether the thermodynamics properties (equation of sates) of 
the matter stay the same as those of the real thermal equilibrium, or not. Either prether- 
malization or real thermalization, the emergence of functional relations among macroscopic 
observables indicates the reduction of the other microscopic degrees of freedom irrelevant 
to the description of the physical properties of the matter. This must be due to a kind of 
randomization processes of these irrelevant degrees of freedom (irrelevant variables) where 
the information of the initial condition is lost by random interactions among them. The 
difference between the real thermalization and pre-thermalization is that, in the case of pre- 
thermalization, there still exist some degrees of freedom which still depends on the initial 
condition. 

If such a prethermalization scenario happens, we can expect some universal behavior 
such as the Boltzmann-Gibbs statistics which leads to the thermodynamic relations. As a 
matter of fact, the Boltzmann-Gibbs statistics is not a unique possibility to be consistent 
with thermodynamics. In 1988, Tsallis proposed a new statistical mechanics where entropy 
is non-additive . So far, extensive studies have been developed and the applications of 
Tsallis' entropy extend to biological systems, high-energy physics, and cosmology. Because 
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FIG. 1: Presence of strong correlation among the two variables in optimisation problem. 

of the non-additivity of the entropy, Tsallis statistics is considered to be important for 
describing correlated systems. From our prethermalization scenario, the Tsallis statistics 
may be attained before the real thermalization is achieved and the later relaxation process 
after the Tsallis statistics may be described by the time evolution of q, which is the parameter 
of Tsallis statistics, and the Tsallis statistics is reduced to the Boltzmann-Gibbs statistics 
for q = 1 0,3. 

To make the idea clear on the effect of correlation on the relaxation time of a prether- 
malization mechanism, let us consider a problem of finding a minimum of a function of 
two variables, as shown in Fig. 1. When the function has a narrow long valley, starting 
from some arbitrary initial point, the usual steepest descent method always first goes to the 
closest point of the valley instead of going directly to the real minimum, as indicated by 
the dashed curve. The narrow valley means that there exists a strong correlation between 
the two variables. In other words, if there exist any strong correlations among variables, 
the minimizing path always tries to satisfy the correlation first (fast relaxation stage), then 
follows downstream to the true minimum along the valley generated by the correlation (slow 
relaxation stage). 

The above consideration suggests that a prethermal state may correspond to a kind of 
narrow valley of the free energy generated by dynamical correlations. If such a local minimum 
stays for a very long time, we may consider the system as in a kind of equilibrium, but the 
microscopic occupation probability is not that of the Boltzmann-Gibbs statistics, due to the 
constraint created by the correlation. Then what characterizes such a local equilibrium? 
We show that there exists a class of dynamical correlation for which the minimisation of 
the energy under the constraints lead exactly to the general form of Tsallis' statistics. In 
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this image, the dynamics of the slow relaxation stage may be described by the evolution 
of q to arrive at the real thermal equilibrium described by the Boltzmann-Gibbs statistics 
corresponds to q = 1. 

Note that in the above analogy, the correlation is rather static one. However the dy- 
namical constraints are not necessarily determined only from the Hamiltonian itself, but 
may depend on the state of the system in a dynamical way. In other words, for a given 
system, depending on the initial condition, the correlation and hence the constraints among 
microscopic degrees of freedom may change. In this sense, it is possible that the value of q 
is not a quantity determined from the beginning even when the system is specified. 

Of course, prethermalization involves a broad range of phenomena, and has been studied 
for a long time in several works |9j . In particular, the relaxation to a quasi-stationary state 
under long-range interactions is strongly suggested in the Hamiltonian mean field model 
[lo| . Here we develop a general phenomenological argument of how the states described by 
the Tsallis nonextensive statistics can be regarded as a prethermalized states. 

II. DYNAMICAL CLUSTERS AND TSALLIS DISTRIBUTION FUNCTION 

In this section, we show that the Tsallis distribution function can be obtained when there 
exists a class of correlations among particles in each single particle energy levels. 

Let us consider a system composed of N particles where the single-particle energy spec- 
trum and the corresponding occupation numbers are given by {e a } and {N a }, respectively. 

Suppose there exists a strong g-body correlation among N a particles in any single-particle 
state a. In general, q is a function of a. However, we assume, for simplicity, the number of 
correlated particles are independent of state a. Due to this correlation among particles, the 
rules to determine the occupation numbers {N a } affects the statistics as described below. 
For a given value of q, the number of possible ways of forming correlated g-body clusters in 
each a state is given by 



for N a ^> 1. Here, we used Stirling's formula for the asymptotic values of T-function. Note 
that q is the average number of particles to form a cluster and does not have to be integer. 
We call these correlated subsystems g-clusters. If the g-clusters are formed with equal a 
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priori probability, the number of such clusters M a should be proportional to the number of 
ways of forming them. Thus 

M a = Const x JV«, (5) 
and the total number of the g-clusters in the whole system is then 



M X>„. (6) 



We may estimate the energy associated with the cluster formation. Since a g-cluster in 
the energy level a carries the single particle energy, q x e a , we have 

E q = qY,M a e a . (7) 

a 

Note that, for q — > 1, the expression above formally leads to the total energy of the system, 
except for a constant factor This fact will be useful for the recovering the usual Boltzmann- 
Gibbs distribution. See the later discussion. 

It is natural to think that, once the number of the correlated particles are fixed by q, the 
most probable configuration will be realized by minimizing the correlation energy, that is, 
the energy associated with these clusters, 

5E q = 0, (8) 

with 

M = Const. (9) 

where the variation should be taken with the occupation numbers, {N a } . This is equivalent 
to the variational problem of the following quantity, 

5j2e«N« = 0, (10) 

a 

with the following constraints, 

N a = N = Const, (11) 

a 

J2K = A = Const. (12) 

a 

That is, we consider the following variational problem, 
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for all JV Q ' s. We get 

(e a + A) N«- 1 - n = 0, 
If we introduce the occupation probability 



then the above solution can be expressed as 



N 



(14) 



(15) 



Pa 



E. 

Xq 



q-1 



(16) 



.1 + s a /X j 

where p! = [i/N q ~ l . Calling A = (q — 1) T, this distribution can equivalently parametrised 
in terms of g-exponential function as, 



p(E) 
p(0) 



1 



9-1 



r q-1 T 

E 

Ex Pq \ ~7p 



(17) 



This form of parametrisation is convenient since it reduces to the the Boltzmann-Gibbs 
distribution function in the limit of q — > 1. 

However, it should be stressed that the thermodynamic properties of the prethermal state 
proposed here are different from that of the Boltzmann-Gibbs statistics. In usual derivation 
of the Boltzmann distribution, we have to maximise the total number of microstates acces- 
sible for a given total energy and number of particles. The variation principle Eq. (1131) in its 
form does not lead to the Boltzmann-Gibbs variational equation for maximising the entropy 
of the system, since for q 1, the second term and the third term degenerates if A and ji 
are finite constant, independent of q. 

On the other hand, we may re-arrange these Lagrange multipliers A and fi in terms of 
new constants, (3 and \x as 



1 



where 



1 



0. 



(18) 



A 



pi = /jL + 



(19) 
(20) 
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and we have used 



The second term 



(21) 



where 



ln 9 (N aj 



a 

i - m- 1 



l-q 



(22) 



(23) 



and in the limit of q — > 1, we have 



In, (JV a )-ln (N a ). 
Thus, for q — > 1, the variational principle, Eq.( Tl8|) leads to 
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£ « N <* ~ g 52 N <* ln Na - A £ ^ 



0. 



(24) 



which is exactly the same equation we get in the Boltzmann-Gibbs statistics when we max- 
imise the entropy 

S BG = k\nW J (25) 

with 



ln W = In 



ln + const. 



(26) 
(27) 



after using Stirling's formula. 

All the discussion above can be done in terms of occupation probability, p a from the 
beginning. The variational problem, Eq. ffTOj) together with the constraints, Eqs. ffTTj) and 
(|T2|) is equivalent to 



with 



52 

a 



(29) 
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and 



where the variation is now taken with respect to {p a } . The corresponding equation to Eq. (118j) 
is 

sl ^e a pl~ -^2p a ln q (p a )+ v^pA = (31) 

\ a a a / 

where the second term is nothing but Tsallis' entropy 

S T sallis = - ^ Pa ln< ? ' ( 32 ) 

a 

Thus, we have shown formally that our system is a physical realization of Tsallis' formulation 
of the generalised statistics. 

Here, there are important remarks. Although we have shown that the energy minimum 
principle of correlated clusters lead to the same results as Tsallis' generalised statistics, the 
physical significance of our example is very different from that of Tsallis' approach. In our 
approach, we imposed the condition to minimise the energy of systems of g-clusters, keeping 
the total number of clusters. This energy is not the total energy of the system, but that of 
a subsystem of the system which are strongly correlated. When we neglect the interaction 
energies, the total energy of the system is the sum of single particle energies so that 

Etotai = '^2 £a ^ a - ( 33 ) 

a 

In this case, if we require conservation of the total energy of the system, we have to add the 
constraint, 

SE total = (34) 

in addition to Eqs. (fTTT) . (|T2l) . By introducing a new Lagrange multipliers, we obtain instead 
of Eq.©, 

p(E) ( 1 - E/e r 



- max 



P (o) \ i + m. 



(35) 



where now p(0) ,T and e max should be determined from the constraints, Eqs.( fTTT) . (fT2l) and 
(133]) . As expected, the energy spectrum has the maximum cut-off energy, £ max which is a 
natural consequence of the conservation of energy. For E s mSuX , the above expression 
reduces to Eq. (fT7|) as expected. 
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In general, the number of the (/-clusters and the coefficient of the probability of the cluster 
formation will depend on a. In this case, we may write the number of clusters as 



M a = K a (36) 

r(<? a + 1) 



The variational problem is then 



5 



q*K a e a p q a a + A ^ K a p q a a -fi^Pa 



a 



0. 



Then the distribution function is 

i 

la -1 



p(E) / K /K E 



p(0) 1 1^^^ 



(37) 



q a -lT 

III. TOY MODEL OF CORRELATED PARTICLES 

To check our prethermalization scenario, let us consider a very simple toy model described 
below. 

1. Initially N particles are distributed as {N , N\, N a , ...} over equally spaced energy 
levels. Let us denote the energy of the i— th particle as Ei. 

2. Choose randomly a pair of particles, say, % and j. 

3. Energies of i and j are updated according to one of the following alternatives: 

(a) The new energies are set to the lower one of Ei and Ej, that is, 

E u Ej -> E[ = E' j = Mm (E h Ej) , (38) 
then choose another k-th (k ^ i,j) particle randomly and attribute 

E' k = E k + Max {E u Ej) - Min{E u Ej) (39) 

to conserve the total energy. 

(b) Change the energies as 

Ei -> E'i = Ei ± AE, 

Ej ^E'j = Ej t AE, (40) 

where AE is the level spacing. Here, if one of E[, Ej becomes negative, then this 
step is skipped. That is, we have to keep always Ei >0. 
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4. The alternatives (a) and (b) are chosen randomly, but the ratio of the average frequency 
of (a) to (6) is kept as a constant, r. 

Here, the collision of type (a) is the process of the formation of a kind of clusters, while 
the collision of type (b) is the energy exchange process and may destroy the formed clusters. 
It is well-known that for r = (no type (a) collision) in the above model, the final single- 
particle spectrum will become the Boltzmann distribution. By the competition of the two 
collision processes, we expect the number of correlated particles will reach some stationary 
value, and, by construction, the energy of a correlated pair always tends to diminish. In 
this way, we may expect that the above system will lead to the situation described by 
Eq. ( Ti~6l) . In Fig. 2, we show the results of simulations, for several values of r from 0.00001 
to 0.1. It is interesting to note that r > leads to a non-Boltzmann distribution which 
is well approximated by the Tsallis distribution, p a = Cj [1 + (q — 1) /^-EJ 1 ^ 1 " 1 ^, where 
(3 = 1/(3 — 2q) (E) and C = (2 — q)/(3 — 2q) (E) are determined from the normalization 
condition and conservation of energy. One parameter fits with respect to q were performed 
and the results are indicated by the continuous curves in this figure. For r — > 0, the spectrum 
converges to the Boltzmann distribution p a — > e~ Ea ^ E ^ as expected. Note that for a one 
dimensional case like the present model, the Tsallis distribution is valid only for q < 3/2, 
otherwise the energy expectation value diverges. For larger values of r, the fitted value of 
q tends to this limiting value, but the distribution begins to deviate substantially from the 
Tsallis distribution in the low energy region. 

It is important to note that these distributions are the stationary and stable ones. We 
confirm that starting from any different initial conditions, the final distribution converges 
uniquely for a given r parameter. Furthermore, the cluster formation through type (a) 
collisions accelerate the speed of the convergence. For example, for r = 0.01 the distribution 
converges 10 times faster than r = case. 

In the above toy model, the time reversal is violated for the collision type (a). However, 
this is not the crucial factor to obtain the non-Boltzmann distribution. We have checked 
this in a more elaborate model which has time-reversal invariance. 
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r= 0.001 
q= 1.25 

100 200 300 r=0 - 00001 
E q=1.00 



FIG. 2: Energy spectrum after a large number of collisions per particle, starting from a distribution 
peaked at (E) = 22.5. The probability of cluster- forming collisions ranges from r = 0.1 (top) to 
r = 0.00001 (bottom). Tsallis distributions fitted to the calculated spectra are also shown, together 
with the corresponding (/-values. 

IV. INTERACTING LINEAR HARMONIC OSCILLATORS 



As an another simple example in which a single particle spectrum deviates from a simple 
exponential form, let us consider N identical harmonic oscillators q(t), described by the 
following Hamiltonian, 

1 / Jn\ 2 1 . 1 

(41) 



2\dt 



\f + \ G fCq. 
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Here 



<?2 (t) 



(42) 



\Qn (t) J 

and C is an N x N symmetric matrix, representing the linear interactions among these 
TV harmonic oscillators. The constant G represents the magnitude of coupling. It is well- 
known that such a system is formally solvable. Let U the orthogonal transformation which 
diagonalize the matrix C, 

UCU T = A, (43) 
where A is the diagonal matrix whose diagonal elements are eigenvalues of C. 



A = diag(X 1 , A 2 , .., Ajv). 



(44) 



Then by introducing the new vector, 



Ut) = Uq(t) 



(45) 



we have 



so that the general solution for £ is given by 



1 



& (t) = — Ai sin (uit) + Bi cos (uit) , i—l,..,N 

Us 



where 



Ui 



(46) 



(47) 



(48) 



and A = (Ai) and B = (Bi) can be determined from the initial condition given for q and 
dq/dt as 

B = Uq(0) (49) 



and 



A = U 



dq(0) 
dt 



(50) 



Eq. (l47p . together with Eqs. (l49l) and (|50|) shows that the whole system never realize the 
thermally equilibrated state, which is a very known fact. On the other hand, a single particle 
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trajectory in a reduced phase space (q, p) is very complicated and even ergodic if the coupling 
matrix C is sufficiently complex. It is then interesting to see how the long-time average of 
single particle phase space distribution of any of qi (t)'s behaves. Let us define 

/ (q, a) AqAq = hm ~ £ £ dt 6 (2) [fa (t) , q t (*)) ; Afi] , (51) 
where is 2— dimensional Heaviside's step function defined by 

e»[(,W,i(.))M = ( n/W1 ' i(,)) ^ (52) 

I otherwise, 

and Af2 is an infinitesimal domain around the phase-space point (q, q) whose volume is 
given by AqAq. In the above, we assume that any of harmonic oscillators are identical, so 
we consider the inclusive single particle distribution of these N oscillators. 

In Fig. [31 we show the results of numerical calculation for the system with N = 100 and 
the symmetric matrix C is created by a homonegeous random number Cy 6 [—1,1]. Also a 
random initial condition, 

% (0) e [0, l] , 
ft(0) e [-1,1], 

has been specified. The result seems an Gaussian distribution in (q,p) plane. To see more 
precisely, we show the average of single particle spectra 



N 

where 



(54) 



and (t) is the single particle energy without interaction term, 



Ei = \v\ + \<£ + \gCuql (55) 



which corresponds to the single particle oscillator energy with renormalised frequency [13| 
In Fig. HI we show the spectrum corresponding to the example of Fig. [3], which in fact, 
the inclusive single particle behavior is that of the thermal equilibrium. In Fig. [H the 
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FIG. 3: Phase Space distribution of single particle states. The abscissa represents the q coordinate, 
and the ordinate represents the q coordinate. 

triangles are those calculated from the expression (I53I) , whereas the dashed line indicate the 
Boltzmann distribution, 

P B (E) oce- E/{E \ (56) 

where (E) = E^ oi jN . We see that the inclusive energy spectrum of harmonic oscillators 
is essentially described by the Boltzmann distribution. However, we should note that if 
we see more carefully individual spectra for distinct oscillators, they fluctuate around the 
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FIG. 4: Energy spectra corresponds to the example of Fig. [3J The numerical results (the triangles) 
are fitted by the Boltzmann distribution function (The dashed line). Thin lines are spectra for 
several individual oscillators. 

Boltzmann distribution (thin lines). 

Such a behavior can be seen in various initial conditions and C matrix. The above 
example corresponds to the situation where all the oscillators couple each other, a typical 
case of long range interaction. We may consider more local interaction, that is, the matrix 
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is semi-diagonal. Let us consider, for example, 

Cij = [R] e~ a ^~ j ^ (57) 

where [R] is a uniform random number between [—1,1] and a is a parameter related to 
the locality of the interaction. In Fig. [5l we show the inclusive spectrum for the case 
of a = 0.2. Surprisingly, the inclusive spectrum (triangles) deviates appreciably from the 
expected Boltzmann distribution. Instead, the form of the spectrum can be fitted very well 
in terms of Exp q (-E/T) function, with q = 1.05 and T = 0.9 (E). 

The reason why the inclusive spectrum deviates from the Boltzmann to Tsallis (Exp q ) 
distribution is not obvious. However, we may note several key aspects indicating this result. 
First, the semi-diagonal nature of the matrix C means that there are more correlations 
among oscillators q than they are uniformly coupled among them. Second, the individual 
spectra for different oscillators fluctuate much more compared to the case of uniformly 
interacting case (Fig. H]). Third, the effective temperature T is smaller than the Boltzmann 
case. The last, but not least, important point is that the inclusive spectrum we are looking 
in this coupled harmonic oscillators does not account for the interaction energy. That is, 
the sum of single particle energies defined in Eq. (l55|) does not give the total energy of the 
system. The cross terms as qiCijqj for % ^ j are not counted. These terms are omitted in 
the definition of single particle energies, expecting that these are "irrelevant" information 
and cancels out for the behavior of single particle energy in average. However, these effects 
become relevant in the presence of correlations. 

Of course, there are many cases where the inclusive spectrum is neither Boltzmann nor 
Tsallis. For example, if we generate as homogeneous random numbers between [0, 1] 
(that is, all CV,- > ) or between [—1, 0] (all < ), the inclusive spectrum are completely 
different. 

It is then interesting to study the variation of the inclusive spectrum with respect to 
the form of the coupling matrix C. However, it is not a trivial task to choose the coupling 
matrix C in an arbitrary manner, since we have to have all the eigenvalues real and positive 
to guarantee the stability of the system. From this reason, we may generate the interaction 
matrix starting from its eigenvalues. First we choose the set of N positive eigenvalues. Then 
we construct a random unitary matrix U to define 

C = U T AU (58) 
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0.0001 1 — ' — 1 — ' — 1 — 1 — 1 — ' — 1 — 1 — 1 

2 4 6 8 10 

E/<E> 



FIG. 5: The energy spectra in the case of the semi-diagonal random coupling. The triangles 
correspond to the calculated inclusive spectrum and dashed-dot line is the fit by Exp g function. 
Thin lines are spectra for several individual oscillators. 

The generation of the random unitary matrix can be done using the Schmidt orthonormal- 
ization method from a set of arbitrary TV linearly independent vectors, {e*;} . 

In Fig. [HI we show the sequence of inclusive spectra varying the coupling constant G 
for the eigenvalues of C randomly chosen between [0,1]. The initial condition are taken 
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FIG. 6: the sequence of inclusive spectra varying the coupling constant G for the eigenvalues of C 
randomly chosen between [0, 1] . 

randomly generated as in the examples before. We can see that all the inclusive spectra are 
consistent to the Boltzmann distribution up to the energy plotted. (For larger energy, the 
spectra deviate from the Boltzmann due to the finiteness of iV (= 100 ). 

On the other hand, in Fig. 0, we show the sequence of inclusive spectra when the 
eigenvalues are very inhomogeneous. In this example, we took eigenvalues as a quickly 
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increasing function of its index, 



A similar behavior of the inclusive spectrum appears when we take 




(59) 



(60) 



where ( is a constant. We found that all of these spectra are very well fitted by the two 
parameter Tsallis distribution, Exp q (-E/T). The negative curvature spectra also appear, 
but they can also be fitted by q < 1. 

It is interesting to observe that other simple 2-parameter function, such as exponential 
of polynomial in E, 



does not work at all. 

Another important fact is that, as we mentioned in the introduction, the correlation 
among the internal degrees of freedom may depend on the initial condition so that the 
value of q can also depend on the initial condition. In fact, the form of inclusive spectrum 
when we change the initial condition changes but still can be expressed by Exp q (—E/T), 
with different q and T. Therefore, the parameter q is rather a kind of state parameter like 
other thermodynamic quantities and not the parameter specified by the Hamiltonian of the 
system. 

V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we discussed the concept of the prethermalized state as the presence of 
strong correlation in the phase space of the system. We argued that those states described 
by Tsallis statistics belong to this category. We show that a very specific class of correlation 
leads naturally to the Tsallis distribution for the occupation probability in energy levels. 

We also show an interesting simple example where Tsallis distribution emerges. We 
studied the inclusive energy spectrum of N coupled harmonic oscillators. Depending on the 
coupling, the inclusive spectrum exhibit clearly the Tsallis distribution. The discussion here 
have been restricted in the framework of classical mechanics. However, the mathematical 
structure of general coupling matrix is equivalent to a quantum mechanical problem to find 



exp (-aE + bE 2 ) 



(61) 
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FIG. 7: The sequence of inclusive spectra when the eigenvalues are very inhomogeneous. 

the eigen frequencies if the vector q is considered as a state vector in an appropriate base. 
The single particle spectra are thus somewhat related to the reduced density matrix defined 
for each individual oscillators. In this sense, the question treated here can also be related 
to the problem of decoherence and locality. 

The emergence of Tsallis distribution is not deterministic, but probabilistic in the sense 
that it depends on the choice of the random unitary transformation. There are specific cases 
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that the inclusive spectrum does not behave neither Boltzmann nor Tsallis, at all. However, 
when we choose the unitary transformation randomly, then the probability of emerging 
Boltzmann, or Tsallis type distribution seems to be dominant. For example, if we go from a 
system where the coupling strength is small and the interaction matrix is uniformly random 
to the system where there exist stronger interactions or correlations, the inclusive spectrum 
starts to deviate from the Boltzmann to Tsallis, statistically. This observation deserves a 
further systematic investigation. It will be interesting if we can formulate the problem from 
the point of variational problem as discussed in this paper. 

We thank fruitful discussion with Carlos Aguiar and Gabriel Denicol. 



[1] C. Bustamante, J. Liphardt and F. Ritort, Physics Today, July (2005) 43. 

[2] See, for example, Hama Y, Kodama T and Socolowski Jr O, 2005 Braz.J.Phys. 35:24-51; P. 

Huovinen and P.V. Ruuskanen, 2006 Ann. Rev. Nucl. Part. Sci. 56, 163; Jean- Yves Ollitrault, 

2008 Euro. J. Phys. 29, 275 and references therein. 
[3] P. Romatschke and M. Strickland, Phys. Rev. D68, 036004 (2003); P. Arnold and J. Lenaghan, 

Phys. Rev. D70, 114007 (2004); M. Strickland, Braz. J. Phys. 37, 762 (2007) and references 

there in. 

[4] J. Berges, S. Borsanyi and C. Wetterich, Phys. Rev. Lett. 93, 142002 (2004). 
[5] C. Tsallis, J. Stat. Phys. 52, 479 (1988). 

[6] See for example, S. Abe and Y. Okamoto (Editors), Nonextensive Statistical Mechanics and 
its Applications, Lecture Notes in Physics, Vol. 560 (Springer- Verlag, Berlin 2001); S. Salinas 
and C. Tsallis (Editors), Braz. J. Phys. 29, No. 1 (1999). Also for the recent development, see 
C. Tsallis, this volume. 

[7] T. Kodama, H.-T. Elze, C. E. Aguiar and T. Koide, Europhys. Lett. 70, 439 (2005). 

[8] H.-T. Elza and T. Kodama, Phys. Lett. A335, 363 (2005). 

[9] See for example, D Jou, J Casas- Vazquez and G Lebon, Rep. Prog. Phys. 62, 1035 (1999), L. 
Berthier, J.-L. Barrat and J. Kurchan, Phys. Rev. E61, 5464 (2000), Th. M. Nieuwenhuizen, 
Phys. Rev. E61, 267 (2000); T. Dauxois, S. Ruffo, E. Arimondo, M. Wilkens (Eds.), Dynamics 
and Thermodynamics of Systems with Long Range Interactions, Lecture Notes in Physics, 
Springer, 2002. 



23 



[10] H. Morita and K. Kaneko, Phys. Rev. Lett. 94, 087203 (2005); F. Baldovin and E. Orlandini, 
Phys. Rev. Lett. 96, 240602 (2006); PH. Chavanis, Euro. Phys. J B 53, 487 (2006) and 
references therein. 

[11] In this paper, we use the natural unit, c = h = 1. 

[12] More generally, Eq.([3|) is valid if the four-energy momentum tensor is traceless (conformal 
theory) and isotropic in the Landau local rest frame. For a hydrodynamic description, we only 
need the thermodynamic relations (equation of state) as Eq.Q. 

[13] Another possible way to define the single particle energy is to take the part of Hamiltonian 
depending on z-th coordinate, E{ = ^pf + |g? + ^gJ2f=i Cijqiqj such that J2iL\ Ei = Et t- 
However, this energy is not positive definite depending on the interaction. 



24 



